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ABSTRACT 

^ ■ We determine the relative rates of short GRBs in cluster and field early-type 

CN ' galaxies as a function of the age probability distribution of their progenitors, 

• Pij) oc r". This analysis takes advantage of the difference in the growth of 

pL| ' stellar mass in clusters and in the field, which arises from the combined effects 

lO ' of the galaxy stellar mass function, the early-type fraction, and the dependence 

of star formation history on mass and environment. This approach complements 
CN ■ the use of the early- to late-type host galaxy ratio, with the added benefit that the 

\Q • star formation histories of early-type galaxies are simpler than those of late-type 

^ ■ galaxies, and any systematic differences between progenitors in early- and late- 

Q\ • type galaxies are removed. We find that the ratio varies from -RciusterZ-RfieW ~ 0.5 

^ ' for n = —2 to ~ 3 for 72 = 2. Current observations indicate a ratio of about 2, 

O ■ corresponding to n ~ — 1. This is similar to the value inferred from the ratio of 

^ ■ short GRBs in early- and late- type hosts, but it differs from the value of n ~ —1 

for NS binaries in the Milky Way. We stress that this general approach can be 
2 ■ easily modified with improved knowledge of the effects of environment and mass 

^ ■ on the build-up of stellar mass, as well as the effect of globular clusters on the 

^ ■ short GRB rate. It can also be used to assess the age distribution of Type la 

•'-j ■ supernova progenitors. 

5r ! Subject headings: gamma-rays:bursts — galaxiesxlusters — galaxies:formation 



1. Introduction 



G amma-ray bursts (GRBs) are divided into two broad classes of short/hard and long/soft 
bursts (IKouveliotou et al.lll993l ). which appear to have different progenitor populations. Ob- 
servations of Type Ic supernovae in association with long GRBs provide a direct confirmation 
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that they arise from the death of massive stars (iHjorth et aLll2003l : IStanek et aLll2003l ). Short 
GRBs, on the other hand, have long been suspecte d to arise frorn. the merger of neutron star 



and/or black hole binaries (NS-NS, NS-BH; e.g . , lEichler et all Il989l : iNaravan et all Il992 



Rosswog fc Ramirez- RuizI I2OO2I : lAloy et al.ll2005l ). During the last year, follow-up observa- 
tions of Swift and HETE-2 short GRBs have provided initial confirmation to this idea, based 



in particular on the l ocalization of some s hort GRBs to elliptical galaxies (IBerger et al.ll2005 



Gehrels et al. 



2005 



Bloom et al.l 2006a), their lack of ass ociation with bright supernovae 



JFox et allboosl : iHiorth et al.lboosi : ISoderberg et"aDl2006h. and their lower energy release 



and wider beaming angles compared to long GRBs (IBurrows et al.l l2006l : ISoderberg et al. 
2006h . 



Despite this observational progress we are still missing a clear understanding of the 
progenitor population, due to the lack of direct observations (circumburst chemical abun- 
danc e s, gravitationa l waves, or a sub-relativistic, radioactive component: iLi fc Paczyhski 
19981 : lKulkarnill2005l ). Thus, statistical studies of the bu rst properties can be highly effec 



five in understanding their progenitor population (s) (e.g.. iGal-Yam et al.ll2005l : iNakar et al. 



2OO5I : iGuetta &: Pirarul2006l ). Recently, IZheng &: Ramirez-Ruia (120061) investigated the ratio 
of short GRBs in early- and late-type galaxies as a constraint on the age distribution of 
the progenitors. Their analysis combines the global star formation rate in each galaxy type 
with a local galaxy stellar mass function, assuming the formation process of short GRB 
progenitors does not depend on any other physical parameters. Since each galaxy type has 
experienced a globally different star formation history, the ratio of bursts in each type is 
predicted to vary as a function of the progenitor age distribution. 

A complementary way to constr ain the progenitor ag e distribution is to use the rates of 
short GRBs in clusters and the field (IBerger et al.ll2006bl : hereafter Paper I). This approach 
takes advantage of the following differences between cluster and field environment. First, the 
galaxy st ellar mass function of clusters is rnore h eavily dominated by massive galaxies than in 
the field ( ICroton et al.ll2005l : iBaldry et al.ll2006l ). Since the star formation history is mainly 
determined by galaxy mass, the overall growth of stellar mass is in turn affected by the large- 
scale environment. Second, the fraction of early-type galaxies is larger in clusters than in 



the field (jDresslerlll980l : IWhitmore et al.lll993l : iBaldry et al.ll2006l ). Finally, there appears to 
be a systematic offset i n the star formation histories of cluster and fi e ld early-type galaxie s 



of ~ 1 — 3 Gyr (e.g., iBernardi et al.l Il998l : iKuntschner et al.l 120021 : iThomas et al.l l2005l ). 



These effects lead to an overall difference in cluster and field star formation histories, which 
combined with the progenitor age distribution, is expected to affect the relative fraction of 
short GRBs in each environment. 



In this paper we quantify these effects and show how the ratio of short GRBs in cluster 
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and field early-type galaxies can be used to understand the age distribution of the progenitors. 
This is part of our on-going systematic study of galaxy clusters hosting short GRBs, using 
multi-slit optical spectroscopy and X-ray observations (Paper I). We find that the current 
observations, albeit with a small number of events, favor P{t) oc r" with ~ — 1. Finally, 
we provide a comparis o n of the systematic effects in this approach and the approach of 



Zheng fc Ramirez-Ruia (120061 ). 



2. The Rate of Short Bursts in Clusters and the Field 



The approach used by IZheng &: Ramirez-Ruia (120061 ) can be generalized to formulate 
the short GRB rate per unit volume in cluster and field early- type galaxies at z ~ 0, using 
the star formation history function, SFH(r) instead of a star formation rate function: 



pt(z=oo) 

Ri = C SFRiir) P{t) dr. 

Jo 



(1) 



Here r is both the look-back time and the time delay of a short GRB progenitor, and i 
designates a cluster or field environment. -P(t) represents the time delay probability dis- 
tribution of the progenitors with a normalization constant C. In the context of NS-NS or 
NS-BH mergers we adopt the standard power-law form, P(t) oc r". We note that NS-NS 
binaries in the Milky Way appear to follow P(t) oc (jChampion et al.l 120041 ). It is thus 
the convolution of P{t) with SFH(r) that determines the relative rate of short GRBs in 
cluster and field early-type galaxies. 

Since the star formation history of early-type galaxies is determined by both the mass 
and environment of a galaxy, the total star formation history of each environment can be 
described in the following manner: 



SFRAr) 



,(M) SFHg,i,i(r,M) dM, 



(2) 



Ml 



where (j)i{M) is the galaxy number density function, SFHgai,i('r, M) is the star formation 
history function of a single galaxy that has a stellar mass, M, and Mu and Mi are the 
appropriate upper and lower mass integration limits (see §3]). 

From SDSS observations, it has been determined that the galaxy number density is well 
described by a double Schechter function (IBaldry et al.ll2006| ): 



(j){M)dM = e-^^/*^*[0*(M/M*)°i + <pl{M/M*Y\ 



dM 



(3) 
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where the parameters logM*, (pl, ai, ^"^^ 0^2 for a cluster environment (logS ~ 1.3) are 
11.06, 0.74, —1.09, 0.07, and —1.5, and for a field environment (log S ~ —0.9) they are 10.44, 
2.7, —0.2, 0.8, and —1.5 (IBaldry et al.ll2006l ). Here, E is the projected density of neighboring 
galaxies, and the specific values are determined from a samp le of ^ 1.5 x 10^ ga laxies in the 
redshift range 0.01 - 0.085 in the SDSS Data Release Four JSaldrv et al.lbooeh . 



We are here only interested in the rates of short GRBs in early-type galaxies, whose 
star formation history is better understood than those of late-type galaxies, and in order to 
avoid any systematic differences between progenitors in late- and early-type host galaxies. 
We therefore need t o modify (/>(M) by the fraction of early-type galaxies in each environment 
tealdrv et aPbooeh : 



/,,,(S,M) = 1 -exp{-[(S/6i)'^ + {M/h,f^]}, 



(4) 
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where the values of the parameters 61, 62? ^3? and 64 are lO'^'^^ Mpc~^ , 0.69, lO^'^'^^ M 
and 0.59. Thus, 0j(M) = 0(M)/r^j(E, M). In a cluster environment, with a higher E 
and systematically larger masses, the early-type fraction is larger than in the field. A plot 
of fr(f){M) for cluster and field environments is shown in Figure [1^. Clearly, the most 
massive galaxies reside preferentially in clusters, while the bulk of the mass in galaxies with 
M < 10^°-5 is in the field. 

Finally, we turn to the star formation history function, SFHgai(T, M). Early-type galax- 
i es in different enyironraents show different s t ar formation histories for a given galaxy mas s 
( iKuntschner et al.l 120021 : iThomas et al.l l2005l : iDe Lucia et al.l l2006l : ISchawinski et al.l |2006| ) . 
Moreover, the star formation history is also determined by the galaxy mass. Following 
Thomas et al.l (120051 ). we use a Gaussian form for the star formation history: 



SFHgai(r,M) 



M 



27rAt 



exp[ 



{r - tpeak)^ 

2(At)2 



(5) 



where the peak of the star formation history function, tp p^ki is determined b y both the 
mass of a galaxy and its environment (Equations 2 and 3 of iThomas et al.ll2005l ). and At is 
determined by the m ass of a galaxy. Here, we assume that the low-density environment of 
Thomas et al.l (120051 ) corresponds to the field, while clusters correspond to the high-density 
environmen10. The overall trend is that less massive galaxies form their stars later (i.e., 
smaller tpeak) and over a wider timescale (i.e., larger At). In addition, tpeak in clusters is 
systematically larger than in the field. These effects are shown in Figure Wp, where we plot 
the star formation histories of galaxies with 10^, 10^, lO^'', and 10^^ Mq in both environments. 



^Following Thomas et al. ( 2005h we adopt the foUowmg cosmological parameters: fi^ — 0.3, f^A — 0.7, 
and Hq = 75 km Mpc"^ 
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The final ingredient is the overall normalization of Equation [31 In its curr ent form this 
equat ion was normalized to a total mass of 10^'' Mq for each environment by iBaldry et al. 



(120061 ) ■ Since we are interested in the overall rate per unit volume, we therefore need to know 
what fraction of th e stellar ma ss is in clusters versus the field at z 'y 0. F rom the study of 
Baldry et al.l (120061 ) . as well as iFukugita et al.l (ll998l ) and lEke et al.l (120(351 ) . it appears that 



about 20% of the total stellar mass is included in cluster environmental, which we adopt 
here. 



3. Constraints on the Age Distribution of Short GRB Progenitors 

To illustrate the combined effect of the trends discussed in the previous section we begin 
by making the simplified assumption that there is a single typical galaxy mass, Mtyp, for 
each environment. This typical mass in turn determines the typical star formation history 
of each environment. In this scenario. Equation [21 can be simplified as: 

SFH,(r) = UMtyp) X SFHgai,i(r, M,,^), (6) 



while Mtyp is determined by (IBaldry et al.ll2006l ): 



log (Mtyp) = 10.73 + O.lSlog E. (7) 

Therefore, Mtyp is ~ 10^°'® Mq and ~ 10^°'^ Mq for field and cluster environments, respec- 
tively. 

We are now in a position to use Equations [H [3l [H [5l and [61 to predict the relative rate 
of short GRBs in cluster and field early-type galaxies. The result as a function of the power 
law index n is shown in Figure [21 We note that the absolute scale in this plot is irrelevant 
since we did not integrate over the full mass function, but it provides insight into the overall 
trend. Namely, as n increases, i.e., as the distribution is more heavily weighted to older 
progenitors, the fraction of short GRBs in clusters increases. This can be understood as the 
combined effect of a systematically earlier star formation episode and a higher typical mass 
in clusters. 

A quantitative determination of the relative rates requires a full integration of Equa- 
tion [21 The results of this integration are shown in Figure [31 Overall, the dependence of the 
ratio of short GRBs in clusters and the field on n is similar to the one found in the simple 



^This corre sponds to the assum ption that cluster environments can be described by logE > 0.5; see 
Appendix A of Baldry et al. ( 2006h for details. 
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case. However, the integration over the full mass function of each environment brings out 
additional trends. For a low mass cutoff smaller than 10^'' Mq, there is sharp downturn in 
the ratio for n < —1. This can be understood from the plots of SFHgai(r) x P(r) shown in 
Figures [T]3 and[T]i. In particular, for n = —2 the short GRB rate in the field is larger than in 
clusters since the field has a higher abundance of low mass (< 10^° Mq) early-type galaxies 
(Figured^), which have young or intermediate stellar populations. For such low values of n 
the ratio is sensit ive primarily to recent star for mation. We note that the same effect is seen 



in the analysis of IZheng fc Ramirez- Ruia (120061 ) for the ratio of early- to late-type hosts. 



In this context {n < —1) the lower mass cutoff in Equation [2] plays an important role. 
In Figure [3] we show the effect of setting Mi = 10^, 10^, and 10^° Mq. As noted above, for 
Ml = 10^° Mq we do not see an obvious downturn because galaxies above this limit do not 
exhibit any obvious recent star formation activity. For the lower values of M;, we find that 
the largest downturn is for Mi = 10^ Mq. The reason for this is evident in Figured^, which 
shows that the largest difference between the field and cluster mass functions is at M > 10^ 
Mq. At lower masses the mass functions converge, leading to a ratio that is ~ 1. 

For n > —1, on the other hand, the rate is dominated by the oldest, and hence most 
massive galaxies (Figure [T]i). In this case the ratio does not depend on Mi, and the predom- 
inance of massive galaxies in clusters, along with the systematically earlier star formation 
episodes, results in an increased fraction of short GRBs in clusters for larger n. In fact, for 
n = 2, we find that there are three times as many short GRBs in cluster early- type galaxies 
as there are in field early-type galaxies. 



Discussion 



We now turn to a comparison of our model with observation of short G RBs. To date 
two short bursts have been localized to clusters at z ~ 0.2: GRB 050509b (iGehrels et al. 



20051 : iBloom et al.l l2006af ) and GRB 050911 (Paperl). We note that in the latter case the 
large error circle prev ents an association w ith a specific cluster galaxy, but the large early- 
type fraction of 80% (IBerger et al.ll2006bl ) suggests that the burst was likely hosted by an 
early-type galaxy. On the oth er hand, only one short GRB has been localized to a field 
early-type galaxy, GRB 050724 jBerger et al.ll2005h . We d o not consider GRB 050813, which 



was hosted by a cluster at a much higher redshift ~ 1.8;lBergeiil2006a[) . and GRB 060502b 



which may be hosted by an early- type galaxy (iBloom et al.l l2006bf ) . but whose large-scale 
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environment has not been fully explored ye^ Thus, the current ratio of short GRBs in 
cluster versus field early- type galaxies is about 2:1, with a large uncertainty due to the 
small number of events. From Figure [31 we find that this ratio cor r espori ds to n ~ — 1. 
This value is lower than n > 3/2 claim ed bv IZheng fc Ramirez- RuizI (120061 ). but is in rough 
agreement with — 1 ^ "^^ ^ found by Berger et al. ( 2006cl ) based on their revised redshift 
distribution with 1/4 — 2/3 of all short GRBs at ^ > 1. 

Clearly, in both methods of estimating the age of short ORB progenitors the uncertainty 
in the inferred value of n is currently dominated by the small number of bursts with a known 
redshift, host galaxy, and large-scale environment type. Since this uncertainty will eventually 
diminish with a larger sample of events, it is interesting to consider systematic uncertainties 
in both theoretical approaches. Our analysis suffers from the somewhat poor definition of 
cluster and field environments. We have used an overall cluster mass fraction of 20%, as 
indicated by several researchers, but this number may range fro m 10 to 30%. Secon d, we 
have used the simplified bimodal star formation history model of iThomas et al.l (120051 ). but 
these authors do not use the same quantitative d efinition of galaxy environment that was 
used for the mass functions by lBaldry et al.l (l2006l ). Since we have used representative mass 
functions, and then scaled the results by the overall mass fraction in clusters and the field, 
this effect should not be significant. Third, the uncertainty in the definition of field galaxy 
environment leads to an overall uncertainty of about 20% in our calculated ratio. Finally, 
since for n < — 1 the ratio depends on Mi, it is essential to understand the appropriate low 
mass limit for early-type galaxies. With the current inferred value of n ~ — 1, however, 
this may not be a relevant issue for short GRBs. 

Similarly, the analysis o f the short GRB rate in early- and late-type galaxies performed 
by IZheng fc Ramirez-Ruia (120061 ) also suffers from systematic effects. First, the ratio is 
affected by the uncertain star formation history at high redshift when early-type galaxies 
formed most of their stars. Second, their estimation of star formation rate does not consider 
environmental effects that appear to be important in current obser vations from SPSS and 
2dF, and which we have accounted for here. Finally, as noted by IZheng fc Ramirez-Ruiz 
(120061 ) ■ it is possible that late-type galaxies have an altogether different age distribution of 
short GRB progenitors than early-type galaxies. This problem is overcome by our method 
since it considers only early-type host galaxies. We note that if both derivations are in fact 
correct, then the estimated values of n can be used to assess any systematic differences of 
short GRB progenitors in early- and late-type galaxies, as suspected to exist for type la 



■^The limit on the X-ray luminosity at the redshift of the putative host galaxy, Lx 6 x 10''^ erg s ^, 
is lower by a factor of eight than that of the cluster hosting GRB 050509b, but it is somewhat higher than 
that of the cluster hosting GRB 050911 (Paper I). 
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supernovae (jSuUivan et al.ll2006l ). 



Future applications of our approach will include the effect of globular clusters, which 
are thought to provide an efficient environment for the production of NS-NS binaries, and 



may account for a s ubstantial fraction of all short GRB progenitors (IGrindlay et al.l 12006 



Hopman et al.ll2006l ). We expect that since the spe cific frequency of globular clusters in- 
creases significantly with galaxy mass (jHarrislll99ll ). an association with globular clusters 
will increase the fraction of short GRBs in galaxy clusters compared to the trend shown in 
Figure O Similarly, our approach can be extended to higher redshift to investigate the evo- 
lution in the fraction of short GRBs in clusters. We expect that the lack of strong evolution 
in the last several Gyr likely makes our analysis applicable out to z ~ 1. However, if some 
short GRBs are in fact associated with clusters at z ~ 2, this presents an opportunity to 
assess any systematic changes in the value of n with redshift. 

We end with the following conclusion. If our current estimate of n > continues to be 
supported by future observations, then this implies that the majority of short GRBs in early- 
type galaxies will occur in clusters. This therefore suggests that short GRBs can provide an 
efficient tool for findin g forming gala xy clusters at high redshift, as already appears to be the 
case for GRB 050813 (lBergerll2006al ). Continued near-IR imaging and optical spectroscopic 
observations of short GRB fields may therefore provide an efficient method for finding the 
highest redshift clusters. 



We thank I. K. Baldry, Y.-T. Lin, M. Brown, A. Dressier, and J. Mulchaey for helpful 
discussions. M.-S. S. acknowledges support from the Observatories of the Carnegie Insti- 
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Fig. 1. — A summary of the model ingredients that determine the relative rate of short 
GRBs in cluster and field early-type galaxies. Panel (a) shows the mass function in each 
environment (solid line: cluster; dashed line: field), including the different early- type frac- 
tions (Equation H]), and an overall cluster mass fraction of 20% (§2]). Panel (b) shows the 
star formation history as a function of environment and galaxy mass (red: 10^^ M©; pur- 
ple: 10^° Mq; green: 10^ Mq; blue: 10^ M©). Bottom panels show the product of the star 
formation history with the short GRB progenitor age distribution function for a power law 
index n = —2 (c), and n = 2 (d). Clearly, a distribution weighted to short merger timescales 
heavily favors lower mass host galaxies (and hence the field), while a distribution weighted 
to long merger timescales favors massive host galaxies (and hence clusters). 
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-2 2 

n 



Fig. 2. — Ratio of the short GRB rate in cluster and field early-type galaxies as a function 
of the age distribution power law index, n, assuming that each environment is described by 
a typical galaxy mass (Q. The scale on the ordinate is arbitrary since we do not consider 
the full mass function, but the overall trend is representative. Since the typical galaxy mass 
is higher in clusters than in the field, the typical star formation epoch is earlier in clusters. 
We therefore expect more short GRBs in cluster early-type galaxies when n is high. 
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Fig. 3. — Ratio of the short GRB rate in cluster and field early-type galaxies as a func- 
tion of the age distribution power law index, n, considering the full mass function in 
each environment. Different colors represent the effect of different mass integration lim- 
its: 10^ < M < 10^2 Mq (black), 10^ < M < lO^^ (blue), and 10^° < M < M© 
(red). The current observed ratio, based on only three events, is about 2, suggesting that 
n ~ - 1. 



